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I.  INTRODUCTION 

The  Martin  Marietta  Corporation,  Orlando  Division,  has 
developed  an  aerodynamic  methodology  and  published  a report  (Ref.  l) 
under  Contract  No.  DAAH01-7^-C-0621  with  the  US  Army.  The  report 
describes  the  development  and  application  of  semiempirical  methods  for 
predicting  aerodynamic  characteristics  of  slender  body  and  body-tail 
configurations . 


In  an  effort  to  evaluate  the  Martin  Methodology,  a computer 
program  was  developed  to  incorporate  the  various  theoretical  and 
empirical  procedures  called  for  in  the  Martin  report.  The  results  of 
this  evaluation  are  covered  in  detail  in  Ref.  2.  As  a follov-up  to 
the  Martin  Methodology  analysis,  a few  modifications  of  the  methodology 
were  accomplished  in  an  attempt  to  correct  certain  disagreements  found 
betveen  measured  experimental  data  and  associated  predictions  derived 
from  the  methodology. 

This  report  is  intended  to  be  a user’s  guide  for  the  use  of  the 
methodology  program.  The  methods  and  procedures  used  in  developing 
the  aerodynamic  methodology  are  covered  extensively  in  Ref.  1 and  only 
those  modifications  to  the  methodology  implemented  by  the  Army  will  be 
discussed  in  this  report. 

II.  PROGRAM  CAPABILITIES 

The  Aerodynamic  Methodology  Program  can  predict  aerodynamic 
characteristics  of  slender  body  and  body-tail  missile  configurations, 
to  include  the  following  predictions: 

1)  Body  alone  normal  force,  center  of  pressure,  and  axial  force. 

2)  Isolated  tail  panel  normal  force  and  center  of  pressure. 

3)  Mutual  interference  effects  that  enable  the  prediction  of 
normal  force  and  center  of  pressure  of  body-tail  combinations. 

The  program  in  its  present  form  was  developed  to  predict  these 
characteristics  for  missiles  with  four  fins  in  cruciform  configuration, 
with  the  fin  trailing  edges  flush  with  the  missile  body  base.  The  fin 
numbering  convention,  as  well  as  the  conventions  for  angle  of  attack 
(a)  and  roll  angle  (<j>),  are  shown  in  Figure  1. 


The  range  of  input  parameters  for  the  program  are: 


1 ) Mach  number 

2)  Angle  of  attack 


3]  Tall  leading  edge 
sweep  angle 


0.8  to  3.0 

0 to  180  for  isolated  components 
(roll  angle  = 0),  0 to  1*5  for 

body-tail  combinations  at  arbitrary 
roll  angles  from  0 to  180  . 

0 to  70  degrees 


II 


3 
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U) 

After-body  length 
(calibers) 

6 to  18 

• 

5) 

Nose  length 
(calibers) 

1.5  to  3.5 

6) 

Tail  taper  ratio 

0 to  1 

7) 

Tail  aspect  ratio 

.5  to  2.0  (aspect  ratio  of  two 
tail  panels  Joined  at  root  chord) 

— r 

III.  DISCUSSION 

Certain  modifications  of  the  original  methodology  were  incor- 
porated into  the  methodology  computer  program.  These  modifications  were 
an  attempt  to  correct  certain  disagreements  between  experimentally 
obtained  data  and  associated  predictions  from  the  methodology  which 
became  apparent  during  an  evaluation  of  the  original  methodology 
(Ref.  2). 

The  modifications  are  limited  to  the  body  alone  normal  force 
prediction  methods  and  involve  the  addition  of  newly  obtained  unpublished 
data.  Specifically,  equation  l6,  page  36  as  depicted  in  Figure  20  of 
Ref.  1,  has  been  eliminated  and  replaced  by  values  of  n (correlation 
factor  for  end  effects)  as  shown  in  Figure  2 of  this  report.  In 
addition,  the  crossflow  drag  coefficient  versus  crossflow  Mach  number 
curve  shown  in  Figure  22a  of  Ref,  1,  has  been  replaced  by  Figure  3 of 
this  report. 


Comparisons  of  predictions  generated  from  the  original  methodology 
and  those  with  the  above  modifications,  compared  with  the  experimental 
data,  are  shown  in  Figures  1+  through  10  of  this  report. 
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Figure  bl.  Modified  methodology. 
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Figure  5a.  Original  methodology. 
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7igure  6b.  Modified  methodology. 
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Figure  7a.  Original  methodology. 
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Figure  8t>.  Modified  methodology. 
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Figure  9a.  Original  methodology. 
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?1  fire  10*.  Original  methodology. 


Appendix  A.  PROGRAM  UTTLI7ATION 


The  proper  format  and  sequence  for  the  Input  cards  are  shown  In  Table 
A-l. 


TART.F  A-l. 


CARP 

SYMBOL  RFFFRF.NCF 

FORMAT (FORTRAN  1 

1 

Title  card  for  each  configuration 

1PA6 

2 

AR , S ,CR , TAMPA , ALE 

SF10, 3 

3 

LNP.IAP.P 

1F10.3 

U 

Al.Fl , AI.FP . PEI AF .MACH , RN, OPTPHT . 
OPTMAC 

!«F10,  .3, FI  3.2,212 

AR 

P 

OR 

I AMP  A 
ALE 
I. NT 
IAP 
P 

ALF1 

ALFP 

PF.LAF 

MACH 

RN 

OPTPHT 


OPTMAC  - 


Tail  aspect  ratio  (two  tail  panels  joined  at  root  chord ) 
Tail  Pemi span  including  body,  (b  + PVP,  inches 
Tail  root  chord  - Inches 
Tail  taper  ratio 

Tail  leading  edge  sweep  angle  - Pegrees 
Nose  length  - Calibers 
After-body  length  - Calibers 
Body  diameter  - Inches 

Beginning  angle  of  attack  for  angle  of  attack  sweep  - Peg. 
Ending  angle  of  attack  for  angle  of  attack  sweep  - Peg. 
Angle  of  attack  increment  for  alpha  sweep  - Pegrees 
Free  stream  Mach  number 
Reynolds  number  - /Ft . 

1 to  calculate  aerodynamic  coefficients  as  a function  of 
roll  angle.  Roll  angle  is  automatically  varied  between  0* 
and  op®  in  increments  of  10*.  No  additional  data  cards 
are  required.  See  Cample  Configuration  P Input  Pats. 

0 to  delete  calculation  of  aerodynamic  coefficients  as 
a function  of  roll  angle 

1 allows  for  multiple  Mach  numbers  for  the  configuration 
and  angle  of  attack  range  determined  by  the  four  cards 
Table  A-l.  Additional  Mach  numbers  are  added  after  the 
four  cards  of  Table  A-l,  one  Mach  number  per  card  using 
FI 0.1  Format 

0 for  single  Mach  number  only 
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MULTIPLE  CASE  CAPABILITY: 


Additional  configurations  can  be  input  simply  by  repeating  the  card 
sequence  of  Table  A-l. 

If  OPTMAO  = 1,  however,  a blank  card  must  follow  the  last  Mach  number 
card  before  the  program  will  accept  additional  configurations. 

Table  A-2  shows  a sample  input  card  arrangement.  Using  this  input 
card  arrangement  the  program  will  predict  aerodynamic  coefficients  for 
sample  configuration  No.  1 for  Mach  No.  equal  .8  and  for  angles  of 
attack  from  0.0  to  i*5.  degrees  in  5 degree  increments.  Since  OPTPHI 
is  0 for  sample  configuration  No.  1,  the  coefficients  will  be  given 
versus  angle  of  attack  only  (roll  angle  equals  0 degrees).  Now, 

OPTMAC  equals  1,  therefore  additional  predictions  will  be  made  for  the 
same  configuration  and  angle  of  attack  range  for  Mach  numbers:  . Q,  1.0, 
1.1,  1.2,  and  1.3.  A blank  card  follows  the  Mach  number  cards  for 
sample  configuration  No.  1,  which  signals  the  program  to  read  additional 
configurations . 

The  program  then  reads  the  cards  for  sample  configuration  No.  2, 
and  computes  the  predictions  of  the  aerodynamic  coefficients  and 
centers  of  pressure  for  this  configuration.  Now,  since  OPTPHI  equals 
1 in  this  case,  these  predictions  will  be  given  as  a function  of  roll 
angle,  for  the  angle  of  attack  range  0 to  Uo  degrees.  Also,  since 
OPTMAO  equals  0 for  this  case,  the  predictions  will  be  made  for  Mach 
No.  equal  2.0  only. 

The  resulting  output  for  these  sample  input  cards  are  found  in 
Table  A-3. 

As  a guide  for  time  and  storage  requirements,  the  following 
figures  for  the  CPC  6600  are  provided: 

Compile  Time  - 1**  seconds 

Execution  Time  - 2 seconds 

Storage  - Less  than  100K. 
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SIMPLE  CONFIGURATION  1 OPTPHIM.  OPTNAC: 
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Appendix  B.  DESCRIPTION  OF  SUBROUTINES  AND  PROGRAM  LISTING 

Table  B-l  lists  the  subroutines  and  gives  their  basic  function. 

A detailed  description  of  each  subroutine  will  not  be  included  in 
this  report.  The  logic  and  procedures  required  of  each  included  method 
are  discussed  in  detail  in  Reference  1.  The  applicable  section  of 
Reference  1 for  each  of  the  Methodology  Program  Subroutines  is  listed 
in  Table  B-l. 

TABLE  B-l.  LIST  OF  PROGRAM  ELEMENTS 

:! 


SUBPROGRAM  NAME 


REF.  1 SECTION 


FUNCTION  OF  SUBPROGRAM 


Controls  input  and  output  to  the  pro- 
gram. Uses  various  subroutines  to 
calculate  aerodynamic  coefficients  and 
centers  of  pressure  for  total  missile 
and  isolated  components. 


Calculates  tail  normal  force  coeffic- 
ient in  presence  of  a cylindrical 
body  as  a function'  of  roll  angle. 
(Angle  of  attack  limited  to  0 to  L 5 
degrees ) . 


Calculates  tail  normal  force  coeffic- 
ient of  an  isolated  tail  (angle  of 
attack  0 to  100  degrees). 


Contains  data  required  by  several 
other  Methodology  subroutines.  Data 
is  taken  from  British  Data  Sheets. 
(Reference  3). 


Calculates  body  alone  normal  force 
coefficient  (angle  of  attack  0 to  180 
degrees ) . 


Calculates  body  alone  center  of  pres- 
sure measured  in  calibers  from  the 
nose.  (Angle  of  attack  0 to  180 
degrees ) . 


Calculates  body  normal  carry  over 
force  coefficent. 


Calculates  tail-to-body  carry  over 
normal  force  center  of  pressure. 


Calculates  the  tail  chordvise  center 
of  pressure  (upper  and  lower  bound)/ 
tail  root  chord,  measured  aft  from 
root  chord  leading  edge. 


SUBPROGRAM  NAME 

REF.  1 SECTION 

FUNCTION  OF  SUBPROGRAM 

YCPTS 

SEC. 5. 2. 3 

Calculates  tail  spanwise  center  of 
pressure  for  isolated  tail. 

XCPTT 

SEC. 5.1. 5 

Calculates  chordwise  center  of  pres- 
sure for  isolated  tail  for  roll  angle 
= 0. 

PHCC 

Required  by  several  other  Methodology 
Subroutines . 

AXIAL 

SEC. 5. 1.3 

Calculates  the  axial  force  coefficient 
for  the  missile  body. 

CSF 

SEC. 5. 1.3 

Calculates  the  skin  friction 
coefficient  for  the  missile  body. 

LOOK1 

One-dimensional  table  look-up  sub- 
routine. 

L00K2 

Two-dimensional  table  look-up 
subroutine. 

LOOK3 

Three-dimensional  table  look-up 
subroutine. 

INDEX 

Required  for  subroutines  L00K1,  L00K2, 
and  L00K3. 

A complete  program  listing  follows. 
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